Thermalization after inflation and production of massive stable particles 
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We discuss thermalization through perturbative inflaton decay at the end of inflation. We find that 
a thermal plasma should form well before all inflatons have decayed, unless all gauge symmetries are 
badly broken during that epoch. However, before they thermalize, the very energetic inflaton decay 
products can contribute to the production of massive stable particles, either through collisions with 
the thermal plasma, or through collisions with each other. If such reactions exist, the same massive 
particles can also be produced directly in inflaton decay, once higher-order processes are included. 
We show that these new, non-thermal production mechanisms often significantly strengthen con- 
straints on the parameters of models containing massive stable particles; for example, stable charged 
particles with mass below the inflaton mass seem to be essentially excluded. 



PACS numbers:90.80.Cq, ll.30.Pb 



TUM-HEP-461/02 



I. INTRODUCTION 

According to inflationary models which were first 
considered to address the flatness, isotropy, and (depend- 
ing on the particle physics model of the early Universe) 
monopole problems of the hot Big Bang model, the Uni- 
verse has undergone several stages during its evolution. 
During inflation, the energy density of the Universe is 
dominated by the potential energy of the inflaton and 
the Universe experiences a period of supcrluminal expan- 
sion. After inflation, coherent oscillations of the inflaton 
dominate the energy density of the Universe. At some 
later time these coherent oscillations decay to the fields 
to which they are coupled, and their energy density is 
transferred to rclativistic particles; this reheating stage 
results in a radiation-dominated Fricdmann-Robertson- 
Walker (FRW) universe. After the inflaton decay prod- 
ucts thermalize, the dynamics of the Universe will be that 
of the hot Big Bang model. 

Until a few years ago, reheating was treated as the 
perturbative, one particle decay of the inflaton with de- 
cay rate Fd (depending on the microphysics), leading to 

the simple estimate Tr ~ (FdMpianck) for the reheat 



temperature |g]Q The reheat temperature should be low 
enough so that the GUT symmetry is not restored and 
the original monopole problem is avoided. In many su- 
persymmetric models there are even stricter bounds on 
the reheat temperature. Gravitinos (the spin— 3/2 su- 
pcrpartners of gravitons) with a mass in the range of 100 
GeV to 1 TeV (as expected for "visible-sector" super- 
particles) decay during or after Big Bang nucleosynthesis. 
They are also produced in a thermal bath, predominantly 
through 2 — !■ 2 scatterings of gauge and gaugino quanta 



@. This results in the boundQ Tr < 10^ - 10^ GeV, in 
order to avoid gravitino overproduction which would de- 
stroy the successful predictions of nucleosynthesis by its 
late decay [Q. 

It has been noticed in recent years that the initial stage 
of inflaton decay might occur through a complicated and 
non-perturbative process called parametric resonance]^, 
leading to an out-of-equilibrium distribution of final 
state particles with energies considerably higher than the 
infiaton mass (the preheating stage) ||9|,[l0|. Parametric 
resonance is particularly efficient for decays to bosonic 
degrees of freedom. For inflaton, 0, coupling to scalars x 
via four-leg interactions K^cji^x^ or three-leg interactions 
hm^4>x^ parametric resonance occurs if the coupling h 
satisfies [fol 



*From now on i\/pianck 
reduced Planck mass. 



= 2.4 X 10^ GeV represents the 



^This bound does not hold in models with gauge-mediated 
supersymmetry breaking, where the gravitino is much lighter 
than sparticles in the visible sector, so that the decay of the 
lightest visible sparticle into the gravitino occurs well before 
nucleosynthesis; nor does it hold in models with anomaly- 
mediated supersymmetry breaking, where the gravitino mass 
exceeds visible-sector sparticle masses by an inverse loop fac- 
tor, so that gravitino decays are sufficiently rapid. We will 
comment on this second scenario later. 

"'Recently, non-thermal production of helicity ±3/2 gravitino 
[^, and helicity ±1/2 gravitino from inflaton oscillations 
have been considered. For models with a single chiral super- 
multiplet, the helicity ±1/2 component of the gravitino is the 
superpartner of the inflaton known as inflatino. The decay 
channels of the inflatino have been discussed in Ref. Q . Also, 
it has been suggested Q , and explicitly shown ^ that in real- 
istic models with several chiral multiplets helicity ±1/2 grav- 
itino production is not a problem, so long as the inflationary 
scale is sufficiently higher than the scale of supersymmetry 
breaking in the hidden sector and the two sectors are coupled 
only gravitationally. 
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h(j) > 



(1) 



where (j) is the amphtude of the oscillations of the infla- 
ton field, and is the frequency of oscillations, which 
is equal to the inflaton mass during oscillations^. The 
explosive decay of the inflaton oscillations via such cou- 
plings comes as a result of two effects, x quanta with 

physical momentum k < ^ hipm^ are produced during 
a very short interval each time that (j) passes through 
the origin. These quanta subsequently result in an even 
faster particle production once their occupation number 
exceeds unity. The combined effect destroys the coher- 
ent oscillations of the inflaton very rapidly and brings 
degrees of freedom to which are coupled to the inffa- 
ton to local equilibrium with it and among themselves 
[ pT[ . After some time the momentum of particles in this 
preheat plasma is redshifted below the inflaton mass and 
another stage of inflaton decay will occur. The nature of 
inflaton couplings determines what exactly happens at 
this stage. It is known that four-leg interactions ft.^0^x^ 
alone do not result in complete decay of the inflaton and 
hence interactions which are linear in the inflaton field 
will be required. If such interactions have a large cou- 
pling, so that is satisfied, the final stage of the infla- 
ton decay will be very quick. This usually results in a 
high reheat temperature which could be problematic in 
supersymmetric models. Moreover, the bulk of the en- 
ergy density may remain in coherent oscillations of the 
inflaton until the final decaying stage even if the inflaton 
coupling satisfies the condition in dU). This happens for 
resonant decay to fermions [ llj where decay products can 
not attain occupation numbers larger than one. It can 
also be the case for bosonic parametric resonance once 
one goes beyond the simplest toy models. For example, 
in the instant preheating scenario |l3| the x quanta may 
not build up a large occupation number, if they quickly 
decay to other fields after each interval of production. 
Also, moderate self-interactions of final state bosons ren- 
ders resonant decay much less efficient [ p"4|Jl5[| . 

It is therefore generally believed that an epoch of (per- 
turbative) reheating from the decay of massive particles 
(or coherent field oscillations, which amounts to the same 
thing) is an essential ingredient of any potentially real- 
istic cosmological model [Q. In what follows we gener- 
ically call the decaying particle the "inflaton", since we 
are (almost) certain that inflatons indeed exist. Note 
also that in a large class of well-motivated models, where 
the inflaton resides in a "hidden sector" of a supergrav- 
ity theory its couplings are suppressed by inverse 
powers of Mpianck, and hence so weak that inflaton de- 
cays are purely perturbative. However, it should be clear 
that our results are equally well applicable to any other 



particle whose (late) decay results in entropy production. 

Inflaton decays can only be described by perturbation 
theory in a trivial vacuum if the density n of particles 
produced from inflaton decay is less than the particle 
density in a thermal plasma with the same energy density 
p, i.e. 



(2) 



Even before all inflatons decay, the decay products 
form a plasma which, upon a very quick thermalization, 
has the instantaneous temperature 



9* ^ -H^rdA^pianck 



1/4 



(3) 



since Fd ~ K^m^. 
to hold only at 7J ~ 



where H is the Hubble parameter and denotes the 
number of relativistic degrees of freedom in the plasmaQ. 
This temperature reaches its maximum Tmax soon after 
the inflaton field starts to oscillate around the minimum 
of its potential, which happens for a Hubble parameter 
Hi < 7710. If infiaton decays at this early time are to 
be described by perturbation theory in a trivial vacuum, 
the constraint (|) should be satisfied already at Hi. This 
implies Tmax < m^, i.e. Fd < m^^/M^unck- Note that in 
chaotic inflation, where (j) ~ -^Pianck initially, the same 
bound follows from the constraint ( ij) 
We note in passing that requiring i 
Fd, where the bulk of inflaton decays take place, leads to 
the considerably weaker constraint Fd < Ti^/Mpianck- 
Indeed, this weaker constraint is satisfied in all poten- 
tially realistic inflationary models we know of that satisfy 
(|l|). If TO^/Mpianck > Td > m3/M2j^^^^, the dynamics 
of the Universe for H will be as in the standard hot 
Big Bang scenario. However, for some period of time with 
H < 7710, the thermalization of inflaton decay products 
would have to reduce the number of particles. This indi- 
cates that in fact one-particle decay may not have been 
the dominant mode for reducing the number of inflatons 
during that period. For the remainder of this article we 
will assume that Tmax is indeed smaller than m^. 

In addition to the thermalized plasma with tempera- 
ture given by eq.(H), there will be inflaton decay prod- 
ucts that haven't been thermalized yet, with energy 
/2[]. During this era the energy density of the Uni- 
verse is still dominated by the (non-relativistic) inflatons 
that haven't decayed yet. The scale factor of the Universe 
a then varies as a oc T~^^^ The Universe remains 

in this phase as long as iJ > Fd. This can have various 



^In this article we assume that is essentially constant 
between ti and t^- 



**We assume that this number remains essentially constant 
for all T>Tr. 

ttHere we consider two body decays of the inflaton. The av- 
erage energy of decay products will be of the same order even 
if the inflaton dominantly decays to three body, or higher, 
flnal states. 
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cosmological implications, for example for Affleck-Dine 
baryogenesis and electroweak baryogenesis [^ . 

Here we discuss the production of massive long-lived or 
stable particles [^-p5|. 

Particles with mass < Tr, and not too small cou- 
pling to the thermal bath of inflaton decay products, 
are in thermal equilibrium and their abundance satisfies 
n ~ at ~ Fd- On the other hand, once all inflatons 
have decayed, stable particles with mass > Tr can 
only be pair-produced by the Wien's tail of the thermal 
spectrum; one might thus naively expect their abundance 
to be suppressed by the Boltzmann factor exp(— Tpj/m^). 
This need not be true if these particles are produced at 
H > Td- Particles with mass Tr < < T^ax and 
not too small coupling have an abundance cx at early 
times. Once the plasma temperature drops below m, 
pair annihilation to and pair creation from lighter species 
keeps the particle abundance at equilibrium so long as 
Tann ^ H. Finally, at T = Tf , the comoving number den- 
sity of the particle freezes at its final value. The physical 
number density of the particle at T = Tr is redshifted 
by a factor of (T^/Tf)^. This suggests that the thermal 
abundance of particles with mass Tr < < T^ax is in 
general only power-law suppressed. This has important 
consequences for the production of heavy (i.e. > Tr) 
stable particles ||2^ in general, e.g. for the Lightest Su- 
persymmetric Particle (LSP) |23| and charged stable par- 
ticles II2J] in models with low reheat temperature. 

To date, almost all studies have considered the creation 
of heavy particles from the scattering of two particles in 
the thermal distribution. On the other hand, particles 
produced from inflaton decay have an energy E ~ 
before they thermalize; as already noted, such very en- 
ergetic particles will exist until inflaton decay completes 
at ~ Fd. We argued above that for perturbative in- 
flaton decay to^ > Tmax- It is thus possible that heavy 
particles are efficiently produced from the scattering of 
such a "hard" particle (with energy E ~ off "soft" 

particles in the thermal bath (with energy T) , or off 
another "hard" particle. Note that this allows to produce 
particles with mass > Tmax- Of course, the hard par- 
ticles eventually come into equilibrium with the thermal 
bath; the competition among different interactions thus 
determines the abundance of heavy particles produced 
through this non-thermal mechanism. Another possibil- 
ity is that the heavy particles are themselves produced 
directly in inflaton decays , in which case their abun- 
dance can be even higher. 

In this article we study these issues. We begin by a re- 
view of thermalization of the perturbative inflaton decay 
products, and verify that a thermal plasma can indeed 
build up well before the completion of inflaton decay, if 
the inflaton decay products couple to some light (or mass- 
less) gauge bosons. We then turn to particle production 
from hard-soft and hard-hard scatterings, and from in- 
flaton decay, estimating the abundance of produced par- 
ticles in each case. It is shown that heavy particles may 
be produced more abundantly than previously thought. 



in particular directly from inflaton decay. We will finally 
close with some concluding remarks. 

II. THERMALIZATION AFTER INFLATION 

After inflation the inflaton field starts to oscillate co- 
herently around the minimum of its potential; at some 
later time it will decay to the fields to which it is cou- 
pled. In the perturbative regime the decay occurs over 
many oscillations of the inflaton field, since Fd <SC m^. 
This means that the oscillating inflaton field behaves 
like non-relativistic matter consisting of a condensate of 
zero-mode bosons with mass m^. The decay rate F^ of 
the oscillating inflaton field is then identical to the total 
decay rate of free on-shell inflaton quanta |^ . Most in- 
flatons have decayed by td — (2/3)Fj^^. By that time the 
bulk of the energy density in the coherent oscillations of 
the inflaton field has thus been transferred to relativis- 
tic particles with an initial energy of order From 
then on the Universe is radiation-dominated, and the en- 
ergy of relativistic particles is redshifted as oc i"^''^, 
where a is the scale factor in the FRW metric. The ther- 
malization of the inflaton decay products then sets the 
stage for the familiar hot Big Bang Universe. 

However, inflaton decay does not suddenly happen at 
td. Rather, it is a prolonged process which starts once 
the inflaton field oscillations commence at i = tj, when 
H = Hi ^ m0. The comoving number density of zero- 
mode inflaton quanta at time i, nc{t), obeys the relation 
nc{t) — ni exp[— Fd(t — ii)], where ni is the inflaton num- 
ber density at tj. For ti < t < the Universe is matter- 
dominated, which implies a (x t^/^^ i.e. — (3/2)t. 
In the time interval between t and t+{2/'i)H^^ = 2t, de- 
cay products with energy ~ /2 are produced and their 
physical number density fih at the end of this interval is 
(assuming ti <C td): 

n;,~ni[exp(-Fdt)-exp(-2Fdt)](^^^ , (4) 

where the last factor comes from the expansion of the 
physical volume. Particles which were produced ear- 
lier have a redshifted energy and number density. The 
spectrum of inflaton decay products has been derived in 
p6| , p7| and it has been shown that, if the inflaton decay 
products do not interact with each other, the number 
density and energy density of the plasma is dominated 
by particles with energy ~ m^/2 in the spectrum. 

Thermalization is a process during which the energy 
density p of a distribution of particles remains constant, 
while their number density n changes in such a way that 
the mean energy E of particles reaches its equilibrium 
value ~ T. For a distribution of relativistic particles 
which consists of ne bosonic and n-p fermionic degrees 
of freedom in thermal equilibrium, and with negligible 
chemical potential, we have p = 7r^/30 (jib + |«f) T'^ 
and n = C{3)/n^ [ub + jnp) Therefore, the 
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ratios pi / E and n^/E are measures for the deviation 
from thermal equilibrium. For perturbative inflaton de- 
cay these ratios are initially less than one, so that ther- 
malization increases the number density and reduces the 
mean energy. Complete thermalization therefore requires 
interactions which change the number of particles to be 
in equilibrium. 

There are three types of interactions which help to 
build up and maintain full (i.e. both kinetic and chem- 
ical) equilibrium: 2^2 scatterings, 2 ^ N scatterings 
{N > 3), and particle decays. Elastic 2^2 scatterings 
redistribute the energy between the scattered particles, 
but play no role in achieving chemical equilibrium. In- 
elastic 2 — > 2 (annihilation) reactions can help to main- 
tain relative chemical equilibrium between different par- 
ticle species, but again leave the total number of particles 
unchanged. For particles with energy E ^ T and number 
density n{E), the rate of 2 ^ 2 reactions with an energy 
exchange AE ~ i? is typically given by F ~ a'^n{E)/ E'^ , 
where a is the relevant coupling constant. On the other 
hand, 2 ^ N scatterings and particle decays increase 
the number of particles and hence play a crucial role in 
reaching full equilibrium. The rate of 2 ^ iV reactions is 
suppressed by additional powers of a. Therefore inelastic 
scatterings with a large squared 4-momentum exchange 
\t\ ~ E^ come to equilibrium later than 2^2 reactions 
do. This would result in a late thermalization and a low 
reheat temperature, if these were the most important 
reactions which increase the total number of particles 
1^,^. One possibility to achieve chemical equilibrium 
more quickly is through the "catal yze d thermalization" 
scenario |26|. As shown in Refs. ||2^,^^, in the absence of 
interactions n{E)/E'^ increases with decreasing E in the 
spectrum of inflaton decay products, due to redshifting 
of "early" inflaton decay products. This suggests that 
a seed of particles with energy E <C m^, which consti- 
tute a tiny fraction of the number density and the energy 
density of the plasma, may thermalize much earlier than 
the bulk of the plasma. The large number of created 
soft particles can then act as targets triggering a rapid 
thermalization of the bulk. However, as pointed out in 
prt , elastic 2 — > 2 scattering of particles in the bulk off 
particles in the seed is efficient; it destroys the seed by 
bringing it into kinetic equilibrium with the bulk. Nev- 
ertheless, catalyzed thermalization can still take place 
if inflaton decay products themselves decay before com- 
ing into kinetic equilibrium with the bulk [ p7[ . However, 
it does not seem very probable that decays of on-shell 
particles alone will be enough for building the chemical 
equilibrium. 

In general, inelastic scatterings which produce rela- 
tively soft particles are the most important processes 
leading to full equilibrium. This has recently been il- 
lustrated in Ref. where the i— channel scattering of 
two matter fermions with energy ~ m^/2 (from infla- 
ton decay) to two fermions, plus one gauge boson with 
typical energy E <^ m^, has been considered. The key 
observation is that (for T = and in a flat space-time) 



t— channel scattering is divergent as AE 0. One thus 
has to choose a physical infrared cut-off in order to esti- 
mate the thermalization rate. A reasonable choice is the 
inverse of the average separation between two particles in 
the plasma, f oc where n is the number density of 

the plasma. This leads to a thermalization rate for hard 
particles (with E ^ m^l2) of the ordeij^ 

Fi„ - n • cr(2 ^ 3) - a^n^/^ (5) 

Once Fin > El^ relatively soft gauge boson with energy 
> T are efficiently created. This increases rt, and hence 
the thermalization rate (H) . Note that the increased tar- 
get density over-compensates the increase of the cut-off 
l^lmin- From then on the number of particles with en- 
ergy ~ T increases faster than exponentially and quickly 
reaches its final value, as pointed out in Ref. [00. In- 
elastic 2^2 reactions then rapidly build up full kinetic 
equilibrium. This suggests that Fjn can in fact be con- 
sidered as thermalization rate, with n being the original 
(pre-thermalization) density, n ~ gi,T'^ / {im^), where g^, 
is the effective number of relativistic degrees of freedom; 
here T is the temperature (^) which the plasma would 
have if it were thermalized. 

This argument assumes that 2^3 scattering can ef- 
ficiently produce "soft" particles with energy > T. The 
production of even softer particles would not help in ther- 
malizing the plasma, since it would not slow down the 
parent particle appreciably, and might even lead to a 
density of such very soft particles which exceeds their 
thermal density. 2 — > 3 reactions with all virtual parti- 
cles (in propagators) having virtualities of order 'n}/^ can 
indeed produce "soft" particles with energy Eg ^ n^/'^, 
if these "soft" particles are nearly coUinear with an in- 
cident hard particle, the emission angle being of order 



*We ignore factors of order log E/T ~ log m^/T where T is 
the temperature of the plasma, which could increase the rate 
by a factor of a few. Note also that the authors of Ref. 
took \t\min ~ , but estimated the target density n from the 
density of particles before thermalization, n ~ T'^/m^. This 
results in a lower estimate of the thermalization rate. 

^The situation is slightly more subtle if the soft gauge boson 
belongs to an abelian gauge group (e.g. photon). As we 
will show shortly, only inelastic scatterings with gauge boson 
exchange in the t— or u— channel results in the estimate in 
eq.(|5|). Such diagrams exist for scattering of a fermion (e.g . 
quark) from a soft non-abelian gauge boson (e.g. gluon) |3l[ . 
On the other hand, in inelastic electron-photon scattering the 
channel diagram has an electron as exchanged particle, and 
hence has a rate much smaller than that of electron-electron 
scattering. However, soft photons with E ^ T annihilate into 
soft e'^e^ pairs at a rate which exceeds (P) by a factor of 
l/a. These soft can then serve as targets for subsequent 
scatterings of hard electrons. The thermalization rate, which 
is set by the rate of the slowest relevant reaction, is then still 
approximately given by (M. 
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yj vp-l'^ I {m^Es). One might wonder whether the emis- 
sion of such a very coUinear particle is in fact physically 
distinguishable from no emission at all. We believe this 
is the case. Note first of all that the lower cut on the vir- 
tualities of all propagators implies that the "coUinear" 
particle is in fact not that close to the emitted one in full 
phase space. Moreover, after a time of order l/Fin there 
will be many such soft particles. Scattering of an almost 
coUinear "soft+hard" pair on a soft particle can yield a 
final state with soft particles being emitted at a large 
angle only if the "soft" particle in the initial state par- 
ticipates; the "hard" particle in the initial state can only 
scatter at a very small angle. In other words, even before 
full thermalization, the plasma allows to physically detect 
the presence of soft, coUinear particles in the "beam" of 
hard particles; the process that "detects" these particles 
also removes them from the "beam" . This happens at a 
time scale significantly shorter than Fin, since one only 
needs 2^2 reactions here. 

In order to check whether the inflaton decay products 
can thermalize before the completion of inflaton decay, 
we compute the maximum temperature T^ax of the ther- 
mal plasma; thermalization occurs before inflaton decay 
completes iff Tmax > Tr. We use eq.(||), with H = Fin 
from d), and Fd = gy^r|/Mpianck, resulting in 

rr rj. i 3 (9*\^l^ Mpianck ^ ^ 

r„.ax-TH U (-j ^J^A - (6) 

\ R / 

which is somewhat highci]^ than the estimate in Ref . : 

rr rr ( 'i 3/4-^^Planck /„x 

V 37710 / 

In any case it is reasonable to expect that the largest 
temperature of the Universe after inflation is between 
the values estimated in Eq.(^) and EqJ^). An interest- 
ing observation is that Tmax, from Eq.(|g), grows esc T^'^ 
Even if is near its COBE-derived upper bound |32] 
of ~ 10^'^ GeV, for a chaotic inflation model, and is 
around 10^ GeV (which saturates the gravitino bound) 
Tniax will exceed Tr if the coupling a'^ > 10~®. This is eas- 
ily accommodated for particles with gauge interactions as 
a will be the gauge fine structure constant. For inflation- 
ary models with smaller m^, early thermalization can be 
realized with even weaker couplings. Recall also that 
we assume perturbative inflaton decays, which requires 
Tmax < Together with eq.(^, taking a < 0.1, 

this implies Ti„ax < 10" (10^) GeV for Tr = 10^ (1) 
GeV; this bound is saturated for = 2Tmax- This 
implies in particular that there will be no "wimpzilla" 



'■The two estimates coincide if Fd saturates its upper bound 
of m^/Af|i,„,k- 



(m^ ^ 10^^ GeV) ||2^ production from thermalized in- 
flaton decay products, even if we use the higher estimate 

i) for Tn^ax. 

For Fin > H >Vd, fresh inflaton decays will keep pro- 
ducing particles with energy E ~ mip/2, but these parti- 
cles quickly thermalize by scattering off a large number 
of soft particles in the thermal plasma. The time scale for 
this process is again given by eq.(|^), where now n^/^ ~ T . 
Inflaton decay effectively completes at ~ Fd when the 
plasma temperature is Tr. From then on the Universe 
is radiation-dominated and its temperature will be red- 
shifted as a~^. 



h h 



h * k 

a) 

_ </'3 

(j)! 




b) 

Fig. 1: Typical scattering diagrams with scalar bo- 
son exchange in the channel for: (a) fermions or an- 
tifermions in the initial and final state, (b) inelastic 2^4 
scattering involving only scalar particles. 

Since the estimate for Tmax crucially depends on the 
rate for inelastic scatterings it is instructive to explore 
this issue a little more deeply. As mentioned earlier, 
for processes with a t-channel singularity, e.g. e~^e~ — *■ 
e~'"e~7 via the t-channel exchange of a photon, one finds 
(Tin oc \t\~^ rather than a-m oc s~^. In fact this is the very 
same singularity which exists in e+e^ e^e~ scatter- 
ing in the i— channel. The appearance of |t|~^ in the 
cross section can be understood in the following way. 
(T receives a factor of from the photon propagator 
while the contribution from the phase space integration 
is (X \t\, thus resulting in cr oc \t\~^ . This implies that 
boson exchange in the t— channel is required in order to 
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have a oc \t\~'^ since the contribution of diagrams with 
fermion propagator will be oc (after phase space inte- 
gration). Moreover, it can be shown that for essentially 
massless external particles the contribution from scalar 
boson exchange is also suppressed. For example, consider 
the t— channel scattering of two fermions via scalar ex- 
change, as shown in Fig. la. A fermion-fermion-scalar 
vertex naturally arises from a Yukawa coupling or, in 
supersymmetric models, from the D— term part of the 
action. Note that scalar interactions flip the chirality of 
the fermion line. However, for massless fermions a flip of 
chirality also implies a flip of hclicity, i.e. a spin flip; this 
is forbidden by angular momentum conservation for for- 
ward scattering, where i ^ 0. As a result, the diagram 
of Fig. la has no t— channel singularity at all.|| 

Finally, it is also possible to have inelastic 2 — > 4 
scatterings, shown in Fig. lb, with only scalar parti- 
cles in external and internal lines. A vertex comprising 
four scalars can for example arise from the Higgs self- 
couphng in the standard model (SM), or from the F— 
and D— term parts of the scalar potential in supersym- 
metric models. However, in this case integration over 
the total allowed phase space only leads to logarithmic 
divergencies in the limit of vanishing external masses. In- 
cluding the squared channel propagator in Fig. lb, the 
integral in question can be written as / dtdM^^dM^Q 1 /P, 
where M^^ = (pi + Pj)'^ is the squared invariant mass of 
the pair z,j = 3,4or5,6. |i| reaches its minimum for 
small M^y-, in which case |t|min — (M^iM^e) / s; the inte- 
grals over and Af|g then only give rise to logarithmic 
singularities, as advertised. Ignoring such logarithmic 
factors, the total cross section will thus again behave as 
a oc s~^. 

The fact that only diagrams with light gauge bosons as 
internal and/or external lines give rise to din oc \t\~^ can 
have interesting implications for thermalization, in sce- 
narios where the gauge group is completely broken in the 
early Universe. This may for example happen in super- 
symmetric models where flat directions in the scalar po- 
tential can acquire a large vacuum expectation value 
(VEV) during inflationfj If gauge bosons have a 

mass nig such that a~^rr^/^< mg < m^, then the most 
important diagrams are those with resonant gauge boson 
exchange in the s— channel, with a-m ^ arrig^. More- 
over, for nig > we will have tiin ~ a'^mT^. This 



^This can also be seen easily by computing the relevant Dirac 
traces. Each fermion line gives rise to a separate trace, which 
is proportional to the scalar product of the 4-momenta of the 
two fermions involved. For massless particles these exactly 
cancel the of the squared scalar propagator. 

**For a perturbative decay, inflaton couplings to other fields 
are so small that the mass induced by the inflaton VEV (if 
any) will be much smaller than the inflaton mass. Therefore 
scalar fields can indeed acquire a large VEV during inflation. 



affects the estimate for the thermalization rate Fin and 
may indeed result in a (much) smaller Tmax- As one con- 
sequence, thermal effects on the flat direction dynamics 
|l9| , ^ may be alleviated. Once flat directions start os- 
cillating their VEV is redshiftcd, due to the expansion of 
the Universe, and the induced mass for the gauge bosons 
rapidly decrease. The flat direction oscillations starts 
when H ^ V" , the second derivative of the scalar poten- 
tial. In models where supersymmetry breaking is commu- 
nicated to the visible sector through gauge interactions 
at a relatively low scale (of order tens of TeV) the poten- 
tial at high field values is exceedingly flat. It is therefore 
conceivable that the flat direction induced masses for the 
gauge bosons, if they completely break the gauge group, 
could delay thermalization until late times. This might 
allow to construct models with low reheat temperature 
even if the inflaton decay width Fd is not very small. 
In the remainder of this article we ignore such possible 
effects and use the (rather generous) estimate in eq.(^) 
for Tiiiax- This allows us to study the most efficient pro- 
duction of massive particles which is possible from the 
plasma of inflaton decay products. 

III. HEAVY PARTICLE PRODUCTION 

We now move on to the issue of heavy particle pro- 
duction. In recent years several mechanisms have been 
put forward for creating very heavy, even superheavy, 
particles in cosmologically interesting abundances. For 
instance particle production could take place in a time- 
varying gravitational background during the phase tran- 
sition from inflationary to the radiation-dominated or 
matter-dominated phase Another possibility is to 

create supermassive particles from preheating p6| . Here 
we will focus on production of very massive particles from 
various processes, including a thermal bath, during per- 
turbative reheating. Note that particle production from 
other sources, if present, would further strengthen the 
bounds which we will derive as they simply add to pro- 
duction from mechanisms discussed here. 

The situation for a weakly interacting massive parti- 
cle (WIMP) X with mass ^ Tr is very well estab- 
lished Such a species is at equilibrium with the 
thermal bath for T > and its number density fol- 
lows (X T^. Once the temperature drops below m^, 
the particle becomes non-relativistic and its equilibrium 
number density becomes n'^ ~ {Tm^)^^^exp{—m^/T). 
Pair annihilation of xs to light particles, occuring at 
a rate approximately given by Fann ~ ce^n^'m~^, pre- 
serves its chemical equilibrium with the bath so long as 
Fann > H. Here is the effective coupling constant of 
X to particles in the plasma. The annihilation rate even- 
tually drops below the Hubble expansion rate, mostly 
due to the fact that n'^ is exponentially suppressed. 
Then the comoving number density of x will be frozen 
at its final value ^ {Tfm^)^^^exp(—m^/T[), where 
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Tf denotes the freeze-out temperature. In the post- 
inflationary era, i.e. for H < Fd, typically Tf ^ (m^/20) 
if a, up to logarithmic corrections |3^. Taking 

a^^ ~ a{m^/Mz)'^ gives a lower bound of a few GeV 
on the mass of heavy, stable neutrinos, known as the 
Lee-Weinberg bound. In the other extreme, the unitar- 
ity bound cr(xx anything) < 47rm~^ provides a firm 
upper bound of about 100 TeV on any stable parti- 
cle that was in thermal equilibrium at any temperature 

The upper bound on might be relaxed consider- 
ably if the initial thermal equilibrium condition for x is 
relaxed, i.e. if Tf > Tr. Even in that case the x particles 
can have been in thermal equilibrium with the plasma 
of SM particles, whose temperature T can significantly 
exceed Tr at sufficiently early times t < l/Fd. However, 
we saw above that at those times the thermal bath did 
not dominate the energy density of the Universe. This 
gives rise to a significant difference from the freeze-out 
situation discussed above, due to the different T depen- 
dence of the redshift factor. In this case x P^^ir annihila- 
tion reactions freeze out at a higher temperature, since 
H ^ v^TV(T2Mpianck) > VS^T^/Mpianck for H > Td. 
Also, is now redshifted oc T*, which is much faster 
than in a radiation-dominated Universe. The situation 
in this case has been investigated in detail in Refs. 22-2^] 
where the relevant Boltzmann equations governing the 
production and annihilation of xs are solved both nu- 
merically and analytically. In Ref. out of equilibrium 
production of x from scatterings in the thermal bath is 
studied and the final result is found to be (the superscript 
"ss" stands for x production from "soft -soft" scattering; 
see below) 



200 

9* 



3/2 



, / 2000TrY 
(x not in equilibrium) 



(8) 



Here fl^ is the x mass density in units of the critical 
density, and h is the Hubble constant in units of lOO 
km/(s-Mpc). We have taken the cross section for x pah 
production or annihilation to be tr ~ a^/m^. Most x 
particles are produced at T ~ m^/A |23j. The density of 
earlier produced particles is strongly red-shifted, while x 
production at later times is suppressed by the Boltzmann 
factor. It is important to note that f2^ is only suppressed 
by [T-^/m^y rather than by exp(— tti^/Tr). We thus see 
that a stable particle with mass 3> lOO TeV might 
act as the Dark Matter in the Universe (i.e. fi^ ~ 0.3), 



if 



2000 Tr • a 



2/7 



Recall that for < 20Tr 



the standard analysis holds as the freeze-out occurs in a 
radiation-dominated Universe. 



* Provided the comoving entropy of the Universe remained 
essentially constant for T <Tf. 



Eq.(||) predicts a relic density that increases with the 
X coupling strength. However, this is true only if the 
X density is always smaller than its equilibrium density, 
which requires ||23| 



1^-1 



1/2 



-^^PlanckT]^ 



(9) 



If this condition is violated, today's x relic density is 
given by 



200 

.9* 



8 • 105 GeV 



(x in equilibrium). 



(10) 



where the exponent a = (1) if XX annihilation pro- 
ceeds from an S {P) wave initial state. The freeze- 
out temperature is now given by X{ = m^/Tf ~ 

\og{0.08g^^^'^alxf^~''MpunckTl/ml), as compared to 

X{ log(0.2(7r"^''^a^a;f ■^""Mpianck/™^) freeze-out oc- 
curs at T < Tr. 

Of course, eqs.(^ and ( p^ ) are only applicable if 
nix ^ 2T,nax, cf. eq.(^) and the subsequent discussion. 
This requires 



Tr 



V200/ 



1/8 M, 



3/8 
Planck 



^ R '"-<* 



< 4a 



V2OO/ ^ Tr 



1/3 



(11) 



in the second step we have used the condition T^ax < 
7710/2. For example, for a — 0.05, Tmax ^ IOOOTr is only 
possible if Tr < 2- 10~-^^Afpianck- Eq-(^ shows that equi- 
librium is more difficult to achieve for larger 771^; this is 
not surprising, since the cross section for xp^^ir produc- 



tion scales like . This means that eq.(g|) will usually 
be applicable in models with large 777^ and comparatively 
smaller Tr , while in the opposite case eq.(^0|) may have to 
be used. We will see shortly that both situations can arise 
in potentially interesting scenarios. Finally, for fixed Tr 
and 777^ , X production from the thermal plasma is maxi- 
mized if the condition {m is saturated. This gives 



j^^s,max/^2 



3 • 10 



25 



200 



9* 



?77x^/pianck 



(12) 



So far we have only discussed x production from the 
scattering of "soft" particles in the thermal distribution. 
However, "hard" particles with energy E ~ ^ T 

are continuously created by (two body) inflaton decay for 
H > Td. These particles eventually thermalize with the 
bath, but this takes a finite amount of time. The pres- 
ence of hard inflaton decay products can therefore affect 
heavy particle production in two ways. Firstly, xs can 
be produced from 2 — > 2 scatterings of a hard particle off 
either soft particles in the thermal bath (if kinematically 
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allowed), or off other hard particles. Moreover, xs might 
be directly produced from inflaton decay. In both cases 
X production may be enhanced and this is the issue to 
which we now turn. 



A. Particle production from hard— soft scatterings 

It is seen from eq.(^) that the energy density in 
inflaton decay products created in the time interval 
[(3/2)iJ~^, 3i7^^], for H > Fd, is comparable to the en- 
ergy density in the existing thermal plasma which has 
a temperature T ^ (iJrdMpianck)^^*- This implies that 
thermalization of these hard particles will increase the 
comoving number density of the thermal bath, roughly 
speaking by a factor of 2, in addition to bringing the 
decay products into kinetic equilibrium. To the accu- 
racy we are working in, we can therefore set the density 
of "hard" inflaton decay products produced during that 
time interval^ to fih ~ p/tu^ ^ O.Sg^T'^/m^. In order to 
estimate the rate of heavy particle production from these 
"hard" inflaton decay products, we also have to know the 
time needed to reduce the energy of these hard particles 
from a value ^ to a value near T. As shown in 

ref. [Q, 2^2 scattering reactions are not very effi- 
cient. The reaction rate is large, but the average energy 
loss per scattering is only ©(T^/m^), giving a slow-down 

time of order [a^T^/m^] ^ (up to logarithmic factors). 
On the other hand, we saw at the end of Section 2 that 
inelastic 2^3 reactions allow large energy losses (in 
nearly coUinear particles) even if all virtual particles only 
have virtuality of order T. The slow-down rate is thus 
again given by eq.(^, where the "target" density is now 
n ~ 0.25,r3:| 



slow 



3a^T 



( — ) 
V200/ 



1/3 



(13) 



Next let us estimate the rate for x pair production 
from hard-soft scatterings. This process is kincmatically 
allowed so long as ET > Arn^, where E is the energy 
of the hard particle so that the square of the center- 
of-mass energy is typically a few times ET. The hard 
particle initially has energy E ~ m^j^ and average (see 
above) number density fih ^ g^.T'^ / {3m^) , just after its 
production from inflaton decay. It looses its energy at 
the rate given in eq.(O). Note that the time required for 



^More exactly, n^, — drih/dt, where only the production 
of hard particles from inflaton decay is included in duh/dt, 
i.e. dufi/dt — r^n^{t). 

^We note in passing that the logarithmic factors cancel in 
eq.(|l^). We included the factor in the cut-ofi' |t|min since 
(almost) all particles in the plasma have non-abelian inter- 
actions with the exchanged particles in the 2 — » 3 scattering 
diagram. 



the hard particle to drop below the kinematical thresh- 
old £',nin — rn^/T depends only logarithmically on this 
threshold. We ignore this logarithmic factor, and simply 
estimate the slow-down time as l/Fsiow On the other 
hand, the rate for x production from hard-soft scatter- 
ings is approximately given by 



(14) 



where we conservatively assume that x's can be produced 
from scatterings of just one species (e.g. electrons). The 
two contributions in ( p^ ) describe 2^2 reactions with 
squared center-of-mass energy ~ m^T and "radiative 
return" 2^3 reactions, respectively; in the latter case 
the hard particle emits a coUinear particle prior to the 
collision, thereby reducing the effective cms energy of the 
collision to a value near m^. This results in the estimateP^ 



^hs 



(T) 



nh 



slow 



V200/ 
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2^6 



(15) 



for xs produced at temperature T . It is clear that be- 
cause of the redshift factor (Tr/T)^ production close 
to Tthr makes the dominant contribution, where Tthr = 
4m^/m0. In order to make a safe (under) estimate we 
choose the temperature Tq = 2rthr for presenting our 
results; note that the x P^i'' production cross section 
at threshold, s = 4to^, is suppressed kincmatically. If 
Tq < Tr, the physical x density at iJ = Fd, i.e. at 
T = Tr, is simply given by eq.(|l5|) with T = 7r. In this 
case the maximal cms energy at Tr is still above 8m^. 
This means that the total x production cross section will 
be dominated by 2 — > 3 "radiative return" reactions, i.e. 
the second, contribution in eq.(15) will usually dom- 
inate. On the other hand, Tq > Tr leads to a physical x 
density at = Fd of order 



10-' (^) 
V200/ 



(16) 



In this case the contribution from 2^3 processes is 
suppressed (by an extra power of a). 

In order to translate the x density at T = Tr into the 
present x relic density, we use the relation (iSl 



m^rixiT^) Tr 



= 6.5 • 10"^ 



T 

now 

200 m^n^(TR) 



9* 



R-^ now 



(17) 



^More precisely, this is the density of x particles produced 
between t and 2t, i.e. during one Hubble time, during which 
T remained approximately constant. 
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where pji is the energy density in radiation, and we have 
used (rii{/i^)now = 4.3-10"^ [III. Our final results for the 
contribution of hard-soft collisions to the x relic density 
are thus: for Tg < 




^ 10" GeV 

(To < Tr) (18) 



where the second term in the last round parentheses 
comes from 2^2 processes. In the opposite situation, 
we have 



200 

9* 



a3 1013 GeV V 
(To > Tr) 



■X J 



where we have again ignored the contribution from 2 — > 3 
processes. 

Two conditions have to be satisfied for our estimates 
(p^)-(pj|) to be apphcable. First, we need Tq < Tmax- 
Note that this constraint is weaker by a factor m^/m^ 
than the analogous constraint for the applicability of 
eqs.(||), (|1^). Second, at temperature T = max(TR,ro)j 
should be sufficiently small that XX annihilation re- 
actions can be ignored. This is true if the xX annihila- 
tion rate is smaller than the expansion rate, a^ny^/m'^ < 
H{T) ~ V5:TV(MpianckT|). If To < Tr, XX annihila- 
tion is thus negligible if 



Tr 



> 



1 ( 9*\ 
3 V200/ 




Planck 



(To < Tr) 
(20) 



Recall that we need > 20rR, since otherwise x's 
would have been in equilibrium at Tr. Moreover, typ- 
ically a^/a ^ 10~^ for weakly interacting particles. The 
condition ( pO| ) will therefore always be comfortably sat- 
isfied in chaotic inflation, where ^ 10~^Afpianck- If 
To > Tr, the condition for xX annihilation to be negligi- 
ble is 



3 V200/ 



T^Mpianck' 



(To > Tr) 



(21) 



Again, this condition can only be violated if <C lO^^^ 
GeV. We thus conclude that in most scenarios with rather 
heavy infiatons the estimates (|l5|)"(^^ are indeed appli- 
cable. 

Finally, unless is quite close to m^, in which case 
To is likely to exceed Tmax, To is well below m^; if 
> 20Tr, To is then also well below the freeze-out 
temperature Tf even if x's once were in thermal equi- 
librium with the plasma. The x density from hard-soft 
scattering can thus simply be added to the contribution 
from soft-soft scattering. 



For a first assessment of the importance of the con- 
tribution from hard-soft scattering, we compare eq.(p^) 
with the maximal contribution (^2|) from soft-soft scat- 
tering. We find that the contribution from hard-soft 
scattering would dominate if m^/m^ < 5 • 10^^ • (a^/a)^ • 
(mj(./TR)^, where we have taken ~ 200 (as in the 
MSSM) . This condition can only be satisfied if ^ Tr 



is not too far below m^, in which case eq.(12) is not 
applicable anyway, since Tmax < "m-x- I^i other words, 
whenever the soft-soft contribution is near its maxi- 
mum, it will dominate over the hard-soft contribution. 
On the other hand, the hard-soft contribution will of- 
ten dominate over the soft-soft one if x never was in 
thermal equilibrium. For To < Tr, this is true if 
a^m^/ra^ < (mx/TR)^(m^/10i6 GeV), which is satis- 
fied unless Tr and are both rather small. Simi- 
larly, for To > Tr, eq.(|l|) will dominate over eq.(§) if 
(19) > 10^3 GeV • a^{700Tji/m^)'^. Recall that the last 



term must be ^0.1, since otherwise the soft-soft contri- 
bution by itself would overdose the Universe. In this case 
the hard-soft contribution will thus dominate in models 
with rather large infiaton mass. Of course, this means 
that cosmological constraints on the model parameters 
Tr and will often be considerably stronger than pre- 
viously thought. We will come back to this point when 
we present some numerical examples. 



B. Particle production from hard— hard scatterings 

If To > Tr , we should also consider x production from 
scattering of two hard particles. Recall that "hard" par- 
ticles are continuously created as long a.s H > T^- Col- 
lisions of these particles with each other can produce x 
pairs if < Note that this constraint is inde- 

pendent of the temperature. On the other hand, once 
a thermal plasma has been established, the density of 
"hard" particles will always be much smaller than the 
density of particles in the plasma. If hard-soft scattering 
at T = Tr can still produce x pairs, it will certainly dom- 
inate over hard-hard production of x's- On the other 
hand, it's also possible that Tq > Tmax, in which case 
hard-soft scattering (and soft-soft scattering) does not 
produce any x particles. 

Note that the rate of x production from hard-hard 
scattering is quadratic in the density of hard particles. 
We can no longer use our earlier approximate solution 
of the Boltzmann equation in terms of the density fih of 
hard particles produced in the time interval from t to 
2t, since the actual density nh{t) at any given time will 
be much smaller than this. Let us first consider the era 
after thermalization, where a plasma with temperature T 
already exists. As noted earlier, a hard particle will then 
only survive for a time ^ l/Fgiow, see eg. (|l3|) . During 
that era, i.e. for Tmax > T > Tr , the production of hard 
particles from infiaton decays and their slow-down will 
be in equilibrium, i.e. the instantaneous density nh{t) — 
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2rrfn0(<)/rsiow, where is the density of inflatons. This 
latter quantity is given by 2rdn^ ~ Hg^,T*/ (Sm^), where 
we have made use of the fact that T doesn't change too 
much over one Hubble time. This leads to a physical x 
density at T = Tr, for x particles produced during one 
Hubble time after thermalization, of order: 



6 . 1027 . f ll.) 

V200/ 



1/2 



(25) 



hh 



a6Mpianckm2 V200/ 



(22) 



We make the following ansatz for the x production cross 
section: 



+ mm 



aa 



X 2 -2/3 
'^x'^ ""plasma 



(23) 



The first term is again the perturbative 2^2 pro- 
duction cross sections, whereas the second term cuts off 
the t— channel propagator at the appropriate power of 
the density of relativistic particles. This second term 
is needed since eq.(p^ shows that x production from 
hard-hard scattering is dominated by early times, i.e. 
high temperatures. In particular, T > is possible, 
in which case no propagator should be allowed to be as 
large as l/m^. 

The fact that the highest temperatures give the biggest 
contribution in eq. (^2|) raises the question what happened 
before thermalization. In this epoch the hard particles 
by definition didn't have time to slow down appreciably 
(except by red-shifting). Moreover, the co-moving infla- 
ton density remained essentially constant during that pe- 
riod. Hence now nh(t) ~ 2rd(t — ti)n0(t), with n^{t) = 

In the last step we have used nip{T^) ~ g*T^J (Sm^) and 
Fd — Y^T^/Mpianck- Introducing = t^n^ (so that 



where cr^j^ is given by eq.(||) with ripiasma ~ 
.g*T^j^x/(3m0). We have checked that for to^ near 10^'^ 
GeV the x density from hard-hard scattering always 
stays sufficiently small for x annihilation reactions to 
be negligible. It is also easy to see that for Tq < 
the contribution ( p5| ) is less than the hard-soft contri- 
bution (|l^). However, the hard-hard contribution will 
exceed the hard-soft contribution (^9|) if Tq > Tr and 
mx^O.l(m0TpMpianck)"^^^; this could indeed be the case 
if is rather close to but well above Tr. 



C. Numerical examples 

As well known, any stable particle must satisfy fi^ < 1, 
since otherwise it would "overdose" the Universe. How- 
ever, in some cases other considerations give stronger 
constraints on the abundance of x- This can happen, for 
example, for unstable massive particles whose lifetime 
is comparable to the age of the Universe ry. Radiative 
decays of such relics which take place after recombina- 
tion (i.e. a.i t > 10^^ s) are tightly constrained by as- 
trophysical bounds on the gamma-ray background [ p9| . 
The tightest bound arises for decays around the present 
epoch |3^: 



(26) 



Another interesting example is that of charged stable 
particles whose abundance is severely constrained from 
searches for exotic isotopes in sea water [Q. The most 
stringent bound on the abundance of such particles with 
is not affected by the Hubble expansion)',' the pro- electriccharge -1 is derived for masses =i 100- 10000 
duction of x particles before thermalization is described 
by 



GeV 



dX^ 
dt 



n^h^ 



< 10" 



-20. 



(27) 



(24) 



where the factor comes from the factor contained 
in nf^{t). For t ^ ti, X^ will thus grow linearly with 
t, which means that the physical density oc 1/t. Of 
course, this behavior persists only until the onset of ther- 
malization, i.e. for t < 1/Fi„, see eq. (|). Using eq.MtJj 
together with the requirement Tmax < it is easy 

to see that the solution to eq.(|2^) at t — 1/Fin always 
exceeds the result (p^). In other words, x production 
through hard-hard scattering is most efficient before ther- 
malization. This is perhaps not so surprising, since dur- 
ing this early epoch, the hard particles have a higher 
physical density (once i 3> ii) and survive longer than 
after thermalization has occurred. Including the redshift 
from Tmax to Tp and using eq.(^^ we arrive at our final 
estimate 



for heavier particles this bound becomes weaker. 

Having mentioned the different cosmological and as- 
trophysical constraints on long-lived or stable massive 
particles, in Fig. 2 we present three numerical exam- 
ples to compare the significance of hard-soft and hard- 
hard scatterings with that of soft-soft scatterings. We 



plot ily/i as a function of my for 



0.01 and 



a = 0.05. The parameters {T^,m^) are chosen as 
(10** GeV, 1013 GeV) (a), (10^ GeV, lO^^ GeV) (b) and 
(3 MeV.lO** GeV) (c), respectively. 

In Fig. 2a, T^ax ^ 2 • 10^° GeV; this explains the 
very rapid fall of soft-soft contributions (which are never 
in equilibrium for the range of shown) for > 
3 • 10^° GeV. We have multipfied the result of eq.(|) 
with exp(— 2m;^/Tniax) ■ exp(2) if m-^ > Tmax, since this 
contribution needs two factors of the "soft" particle den- 
sity. The exact form of this exponential cut-off is de- 
batable, but it is clear that the soft-soft contribution 



10 




should decrease exponentially, and hence quickly become 



irrelevant, once m-y > T„ 



In the present case, the 
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FIG. 2. The relic density x particles would currently 
have if they are absolutely stable is shown as a func- 
tion of TUx for couphngs a = 0.05, ~ 0.01, and a) 
iTn,m^) = (10* GeV, 10^=* GeV), h){10^ GeV, 10^^ GeV), 
c)(3 MeV, 10* GeV). The soft-soft, hard-soft and hard-hard 
contributions are shown by the dotted, short dashed and long 
dashed curves, respectively, while the solid curves show the 
sum of all three contributions. 



soft-soft contribution is overwhelmed by the hard-soft 
one even prior to this sharp decrease (which makes the 
exact form of the cut-off irrelevant). This is because 
To < Tr as long as < I.l • 10^" GeV, in which 
case hard-soft scatterings which occur at the very last 
stage of inflaton decay are the dominant source of x 
production. However, the hard-soft contribution itself 
drops very quickly for > 10^^ GeV, since Tq exceeds 
Tmax in this mass range. We have multiplied eq.(p^) by 
exp(-To/r,nax) • exp(l) if To > Tmax- Once again the 
exact form of this cut-off is not very important, since 
the hard-hard contribution becomes dominant above this 
point until it is kinematically suppressed, and then for- 
bidden, for w 5 • 10^^ GeV. It is interesting to note 
that hard-hard scatterings can very efficiently produce 
"wimpzillas" up to this kinematical cut-off. In fact, for 
this choice of Tr,, and the overclosure bound (for a 
strictly stable x) and the astrophysical bound (for a late 
decaying x) require that must be bigger than 
otherwise hard-hard scatterings will produce x unac- 
ceptably large abundances. Recall, however, eq.(||) and 
eq.(|2^) which show that fi^*^ cx Tr for fixed m^; reducing 
the re-heat temperature by a factor of 2 would there- 
fore result in an acceptable x relic density if Tq > Tmax- 
In contrast, even without the exponential cut-off, the 
soft-soft contribution would have satisfied the overclo- 
sure bound for > 5 • 10^° GeV, i.e. for the parameters 
of Fig. 2a the hard-hard contribution raises the bound 
on by two orders of magnitude. 

The lower reheat temperature in Fig. 2b leads to a 
lower value of Tmax — 6 • 10'' GeV. In this case soft- 
soft production of x pairs actually were in equilibrium 
for < 2 • lO'' GeV, but this contribution again falls 
very sharply once > Tmax, where the hard-soft con- 
tribution becomes dominant. A kink in the hard-soft 
contribution appears at ~ 4 • 10* GeV, where Tq 
starts exceeding TrQ. The curve flattens out just be- 
fore the kink since here the 2 — > 2 contribution to the 
production cross section become important. The hard- 
soft contribution decays rapidly for > 8 • 10^ GeV, 
where Tq > Tmax- Then the hard-hard contribution nat- 
urally takes over and quickly becomes the only source 
until the kinematical cut-off at 5 • 10^^ GeV. The over- 
closure bound now requires > 10* GeV (by a funny 
coincidence the soft-soft and hard-soft contributions are 
comparable around this lower limit). On the other hand, 
for an unstable x the astrophysical bound in cq.(p^) re- 
quires To > Tmax, i-e. > 10^° GeV, if ~ tu (hard- 
soft and hard-hard contributions are again coincidentally 
comparable around this limit) . Notice that the hard-soft 



*This kink is also barely visible in Fig. 2a at 
GeV. 
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contribution have increased this limit by about a factor 
of 10. It is also observed that now the hard-hard con- 
tribution does not suffice to make x interesting Dark 
Matter particle for any value of m^. This again follows 
from the behavior of this contribution. 

Finally, Fig. 2c displays the results for Tr = 3 MeV, 
near the lower limit required for successful nucleosynthe- 
sis; we also chose a smaller inflaton mass = 10^ GeV 
so that this very low reheat temperature will be achieved 
more naturally. In this case T^ax — 600 GeV. A main 
difference from the previous cases is that now the soft- 
soft contribution is very small even for masses as low as 
~ 1 GeV. The reason is that the abundance of x par- 
ticles, produced from soft-soft scatterings, is determined 
by their annihilation which is at equilibrium in this case 
if < 200 GeV. As a result the soft-soft contribution 
is subdominant for almost all values of m^; as usual, it 
cuts off sharply once m^^ > Tmax- Again a kink in the 
curve depicting the hard-soft contribution is recognized 
at ~ 200 GeV, above which Tq > Tr. Even the 
hard-soft contribution essentially vanishes for to^ > 10^ 
GeV, where Tq exceeds Tmax- The hard-hard contribu- 
tion is now completely negligible since Tr, is very small. 
The only meaningful constraint in this case is that from 
eq.(|2^). It is seen that the hard-soft contribution to x 
production again increases the lower bound on by 
about one order of magnitude, to about 10^ GeV, com- 
pared to the bound derived in Ref. ||2^, which only in- 
cludes the soft -soft contribution. Naturalness arguments 
indicate that the mass of the lightest superparticle (LSP) 
should not exceed 1 TeV; this argument is independent 
of whether or not the LSP is charged. The constraint 
of Fig. Ic is not compatible with this bound. Indeed, 
we find that once the hard-soft contribution is included, 
even for Tr = 1 MeV, = 1 TeV is compatible with 
the bound (|^) only if < 5 • 10^ GeV. This indicates 
that scenarios with a charged LSP would require a model 
of low scale inflation. 

One comment is in order before moving on to the next 
subsection. The production of x particles from soft-soft 
scatterings is independent of the main channel through 
which inflaton decays, so long as the decay products ther- 
malize sufficiently rapidly (i.e. such that < Tmax)- 
For example, assume that (for some reason) the inflaton 
mostly decays to leptons while x can only be produced 
from the scattering of quarks. In this case remains es- 
sentially unchanged since electroweak gauge interactions 
bring quarks into equilibrium with a thermal bath cre- 
ated by leptons. On the other hand, il^ and fl^ will be 
suppressed since a hard lepton cannot produce x from 
2^2 scattering off soft particles in the thermal bath 
and/or off another hard lepton. In this scenario x pro- 
duction from hard-soft or hard-hard processes can still 
occur via 2 — s- 4 scattering reactions; if qq XX is al- 



lowed[], so are £q iqxx and £~^£~ qqxX- The corre- 
sponding contributions fl^ and Sl^ will then be smaller 
than our estimates in eqs.(|l^), ( p^ and ( p5| ) by several 
orders of magnitude. This possibility may be relevant to 
the production of exotic electroweak gauge singlet parti- 
cles with large masses, thus relaxing constraints on the 
model parameters in that case. However, the estimates 
for the LSP and charged stable particle examples will re- 
main unaffected since these species can be produced from 
2^2 scatterings of all standard model particles. 



D. Particle production from inflaton decay 

We now discuss the direct production of x particles in 
inflaton decay whose importance has recently been no- 
ticed [Q. Most inflatons decay at T ~ Tr; moreover, 
the density of x particles produced in earlier inflaton de- 
cays will be greatly diluted. Since inflaton decay con- 
serves energy, the density of inflatons can be estimated 
as ~ O.3.g*T^/TO0. Let us denote the branching ratio 
for — > X decays (more accurately, the average number 
of X particles produced in each (j) decay) by B{(j) x)- 
The X density from (p decay can then be estimated as 
PI: 



^docay^2 



2- 10^5(0 ^X) 



Tr 



m0 1 GeV' 



(28) 



Eq.(p8|) holds if the x annihilation rate is smaller than 
the Hubble expansion rate at T ~ Tr , which requires 



M, 



Planck 



m^J V200/ 



(29) 



This condition will be satisfied in chaotic inflation models 
with ^ 10~^Mpianck, if niy- is large enough to avoid 
overclosure from thermal x production alone. It might 
be violated in models with a light inflaton. In that case 
the true x density at Tr can be estimated by equating 
the annihilation rate with the expansion rate: 



Qmax 
X 



5 • 10^ 



ml /200 
al (1 GeV) • AfpianckTR \ ^ 



1/2 



(30) 



This maximal density violates the overclosure constraint 
< 1 badly for the kind of weakly interacting (ofj^^O.l), 
massive {m^ 3> Tr and rrix 1 TeV) particles we are 
interested inj^ For the remainder of this article we will 
therefore estimate the x density from inflaton decay using 
eqM 



*X =JL,if X is a Majorana particle. 

^Eq.(pd) describes the maximal x density if x decouples at 
T r^Tn. It is not applicable to WIMPs decoupling at T < Tr. 
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We now discuss estimates]^ of B{(j) — > x)- This quantity 
is obviously model dependent, so we have to investigate 
several scenarios. The first, important special case is 
where x is the LSP. If is large compared to typical 
visible-sector superparticle masses, (f) will decay into par- 
ticles and superparticles with approximately equal prob- 
abilityg. This can be illustrated for two possible cases: 
when the main inflaton decay mode is via a superpoten- 
tial coupling, and for a gravitationally decaying inflaton. 
In the first case consider the inflaton superfield com- 
prising the inflaton and its superpartncr inflatino 4>, 
and a chiral supermultiplct 5*, which comprises a com- 
plex scalar field ip and its fcrmionic partner -0, with the 
superpotential coupling 

D im^$2 ^ i/i^**'. (31) 

This superpotential generates the following terms 

h^ipipil) ; h^m^i^cjfip'il), (32) 

in the Lagrangian. It is easily verified that </> decays to 
both of ip and "0 (which in turn decay to matter particles 
plus the LSP) at a rate given by F = {h'^/8n)m^. 

Now let us turn to the case of a gravitationally de- 
caying inflaton. As a simple example, consider minimal 
supergravity where the scalar potential is given by | |4^ ] 



V 



f dG dG 

\dip^ dip* 



Planck 



-^^Planck- 



Here G is the Kahler function defined by 



-""Planck 



log 



Planck , 



(33) 



(34) 



where tfi are the scalar fields in the theory. The inflation 
sector superpotential looks like W ~ {l/2)m^{<^ — v)^ 
around the minimum of the potential where v denotes 



■'■In Ref. simultaneous production of x from the thermal 
bath and from direct inflaton decay has been studied and a 
bound on the branching ratio for such a decay mode is derived. 
Our aim here is to show that higher order processes naturally 
provide a decay mode, even if x is not directly coupled to the 
inflaton. 

^The non-thermal production of LSP has also been consid- 
ered in other contexts. For a neutralino LSP these include 
its production from the decay of Q-balls , in cosmic string 
decays as a possible solution to the observational conflicts of 
WIMP cold dark matter Q, in decay of the Polonyi fleld Q, 
and in moduli decay in connection with models of anomaly- 
mediated supersymmetry breaking |Q and models with in- 
termediate scale unification Non-thermal production of 
axino LSP, as Dark Matter, from neutralino decay has been 
considered in 146|. 



the inflaton VEV at the minimum^ The superpoten- 
tial also includes the familiar Yukawa couplings for the 
matter sector, e.g. ft.„H„Qu where Q and u denote the 
superfields containing the doublet of left-handed and the 
singlet of right-handed (s)quarks, respectively, and H^, is 
the superfield which contains the Higgs doublet giving 
mass to the up-type quarks. Then it is easy to see that 
Eq.(^3|) leads to the following term (among others) 



Planck 



(35) 



in the Lagrangian for inflaton decay to three scalars, in- 
cluding two superparticles^ This implies that the rate 
for production of particles and sparticles through infla- 
ton decays into three light scalars is approximately the 
same. 

Once one goes beyond minimal supergravity, the infla- 
ton can also decay to gauge fields and gauginos. Consider 
the case where the gauge superfields have nonminimal ki- 
netic terms, as a result of Afpianck suppressed couplings 
to the inflaton, in the following form 



fap — 



1 + a 



Planck 



+ h 



iV/, 



Planck 



Then to the leading order one finds the term 



Planck 



(36) 



(37) 



for the inflaton coupling to gauge fields, where a repre- 
sents the relevant gauge group index. This results in 
inflaton decay to a pair of gauge quanta at the rate 
F '^^'^^/-^pianck- corresponding term from the 

kinetic energy of the gauginos results in a derivative cou- 
pling and hence a decay rate which is suppressed by a 



factor m|/m^, where rrig is the gaugino mass 44 1. How- 



ever, there exists another term in the Lagrangian |4^ ] 
which is responsible for gaugino mass from supersymme- 
try breaking by the inflaton energy density: 



(38) 



where F,^ is the i^-term associated with the inflaton su- 
perfield, given by ~ m^(j). It is easy to see that the term 
in (pq), to the leading order, results in inflaton decay 



*In the absence of any superpotential coupling to other mul- 
tiplets which provides a linear term in 0, a non-zero v is re- 
quired for the inflaton decay to take place. 

^The main mode for decay to three scalars is to Hu and 
two squarks of the third generation, because of the large top 
Yukawa coupling. 
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to two gauginos at the rate F ~ "^"^^/-^pianck' which 
is coniparable to that for inflaton decay to two gauge 
quanta^. 

Moreover, aU superparticles will quickly decay into the 
LSP and some standard particle(s). As long as > 
Tr, the time scale for these decays will be shorter than 
the superparticle annihilation time scale even if ~ 
0.1. As a resuh, if x is the LSP, then B{(j> ^ x) - 1, 
independently of the nature of the LSP. 

Another possibility is that the inflaton couples to all 
particles with more or less equal strength, e.g. through 
non-renormalizable interactions. In that case one ex- 
pects B{(j) ^ x) ~ 1/5* ~ 1/200. However, even if (j> has 
no direct couplings to x, the rate ( p8| ) can be large. The 
key observation is that x can be produced in (j) decays 
that occur in higher order in perturbation theory when- 
ever X can be produced from annihilation of particles in 
the thermal plasma. In most realistic cases, (f) — > f fxx 
decays will be possible if x has gauge interactions, where 
/ stands for some gauge non-singlet with tree-level cou- 
pling to (j). A diagram contributing to this decay is shown 
in Fig. 3. Note that the part of the diagram describ- 
ing XX production is identical to the diagram describing 
XX ^ ff transitions. This leads to the following esti- 
mate: 



B{cj) ^ x)4 



a 



4a, 



4m^ 



967r3 



2m^ 

TO0 



(39) 



where C4 is a multiplicity (color) factor. The phase space 
factors have been written in a fashion that reproduces the 



correct behavior for 



as well as for 



0. 



This estimate provides a lower bound on i?((/) — » x) un- 
der the conditions assumed for our calculation of fl^ 
and whenever a primary inflaton decay product can 
interact with a particle in the thermal plasma, or with 
another primary decay product, to produce a XX pair, 
— > X four-body decays must exist. It is easy to see 
that the contribution ( P8|) to the x density from infla- 
ton decay will exceed the hard-hard contribution (^5|) , if 
Tmax can be estimated from eq.(^). If Tq > Tr , the de- 
cay contribution will exceed the hard-soft contribution 
if > 0.3m^Tj^/a^, which is true in almost all models 
that avoid overclosure from thermal x production alone. 
Even if Tq < Tr, the decay contribution will dominate 
over the hard-soft contribution if 



> 150 



Tr 



Tr 



(40) 



where the first and second term in the square bracket 
describe x pair production from 2 — > 3 and 2^2 pro- 
cesses, respectively; see eq.(p^. This condition can be 
mildly violated, i.e. in some cases the hard-soft con- 
tribution may exceed the decay contribution. For ex- 
ample, in the scenarios considered in Figs. 2, as long as 
< m0 we find f^^ecay _ 6C4-105tox/(101° GeV) in a); 
- 0.6C4m^/(10^ GeV) in b); and ~ 2 • lO-^^dmJGeV 
in c). For C4 = 1, there is a narrow range of m^- in 
Figs. 2a,b where 17^°^^^^ < ^l^; however, in this range 
fi^ exceeds the upper bound of 1 significantly. We thus 
conclude that the decay contribution to fl^ will usually 
dominate over nonthermal x production from inflaton de- 
cay products if four-body </> — > x decays exist. 




Fig. 3: Sample diagram for x production in four-body 
inflaton decay. 

Occasionally one has to go to even higher order in per- 
turbation theory to produce x particles from (j) decays. 
For example, if x has only strong interactions but (j) only 
couples to SU (3) singlets, xX pairs can only be produced 
in six body final states, (/) ~^ fflQXX- ^ representative 
diagram can be obtained from the one shown in Fig. 3 
by replacing the x lines by quark lines, attaching an ad- 
ditional virtual gluon to one of the quarks which finally 
splits into XX- The branching ratio for such six body 
decays can be estimated as 



2 ^,2 



Bi^^xh 



1.1 • 10^ 



1 - 



4m^ 



1 - 



2m^ 



(41) 



■'•Note that moduli decay to gauginos through the term in 
( ^il) will be suppressed j^. The reason is that the superpo- 
tential is (at most) linear in the moduli superfield, so that the 
corresponding _F— term does not contain the moduli field. In 
contrast, in almost all inflationary models the dependence of 
the superpotential on the inflaton superfield is quadratic or 
higher. 



Another example where xX pairs can only be produced 
in (f) decays into six body final states occurs if the in- 
flaton only couples to fields that are singlets under the 
SM gauge group, e.g. right-handed (s)neutrinos i^r Q. 
These (s)neutrinos can emit a virtual Higgs boson, which 
can split into a top quark-antiquark pair; one of which 
can emit a virtual gluon, which in turn splits into a 
(strongly interacting) xX pair. In this scenario the factor 
in eq. ( ^l| ) would have to be replaced by the combina- 
tion of Yukawa couphngs /i^^/i(/(167r^). If 2m^ < m^j^, 
XX pairs can already be produced in four body final states 
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from vji decay. The effective — > x branching ratio 
would then again be given by eq.(^9|), with replaced 
by m^j^ in the kinematical factors. 

Finally, in supergravity models with explicit (su- 
persymmetric) x mass term there in general exists 
a coupling between and either x itself or, for 
fermionic x, to its scalar superpartner, of the form 
a (m^m^/Mpianck) </'*XX + h.c. in the scalar potentia]||; 
this term is completely analogous to the one shown in 
( ^)P*| . A reasonable estimate for the coupling strength 
is |49[ a ~ I'/Mpianck, unless an i?— symmetry suppresses 
a. Assuming that most inflatons decay into other chan- 
nels, so that Fdccay ~ ^/g*T^ / Mpii^ack remains valid, this 
gives 



Bi^^x) 



167r^/57Mpianck7R 



(42) 



The production of x particles from inflaton decay will 
be important for large and large ratio m^/Tn, but 
tends to become less relevant for large ratio m^/m^. 



Even if 



X production from the thermal 



plasma (H) will be subdominant if 



B{ 



X) /lOOTRym^lTeV 



> 



(43) 



The first factor on the r.h.s. of ( ^ ) must be ^ 10~® in 
order to avoid over-production of x from thermal sources 
alone. Even if x decays only occur in higher orders 
of perturbation theory, the l.h.s. of ( ^ ) will be of order 
10"'* (10~^°) for four (six) body final states, see eqs.(39), 
(^); if two-body tree-level (f) ^ XX decays are allowed, 
the l.h.s. of (^) will usually be bigger than unity. We 
thus see that even for ~ 10^'^ GeV, as in chaotic 
inflation models, and for ~ IO'^Tr, x production from 
decay will dominate if > 10^ (10^°) GeV for four (six) 
body final states; this agrees with the numerical results 
shown in Fig. 2. As a second example, consider LSP 
production in models with very low reheat temperature. 
Naturalness arguments imply that the LSP mass should 
lie within a factor of five or so of 200 GeV. Recall that in 
this case B{(t) — + x) = 1- Taking a^^ ^ 0.01, we see that 
X production from decay will dominate over production 
from the thermal plasma if < 6 • 10^ GeV for Tr = 1 
GeV; this statement will be true for all < 10^^ GeV 
if Tr < 100 MeV. 



^This term also induces an ^-term from supersymmetry 
breaking by the inflaton energy density. 

"Note that a superpotential mass term for the x multiplet 
will be allowed under a Z2 discrete symmetry. For a contin- 
uous symmetry one requires two multiplets xi and ^2 with 
opposite charges. 



In 1^21 showed that the decay contribution (|2^) by 
itself leads to very stringent constraints on models with 
massive stable x particles. In particular, charged stable 
particles with mass below ^ 100 TeV seem to be ex- 
cluded, unless > to^/2. In case of a (neutral) LSP 
with mass around 200 GcV, the overclosure constraint 
implies to^/Tr > 4- lO^'', i.e. a very low reheat tempera- 
ture, unless X was in thermal equilibrium below Tr; recall 
that B{(t) — > x) = 1 in this case. Finally, if ~ 10^'^ 
GeV a "wimpzilla" with mass ~ 10^^ GeV will be 
a good Dark Matter candidate only if it has a very low 
branching ratio, B{<f) ^ x) ^ ^ ■ 10~® GeV/TR, i.e. if its 
couplings to ordinary matter are very small. 

Our calculation is also applicable to entropy-producing 
particle decays that might occur at very late times. If x 
is lighter than this additional </>' particle all our expres- 
sions go through with the obvious replacement — > 0' 
everywhere. More generally our result holds if (j) de- 
cays result in a radiation-dominated era with Tr > m^i . 
If (j)' is sufficiently long-lived, the Universe will eventu- 
ally enter a second matter-dominated epoch, (j)' decays 
then give rise to a second epoch of reheating, leading 
to a radiation-dominated Universe with final reheating 
temperature Tr^, and increasing the entropy by a fac- 
tor m^'/TRj.. This could be incorporated into eq.(28) by 
replacing Tr T^T^Jm^i > Trj. 

In a similar vein, consider the production of a neutral 
LSP from gravitino decay in gravity-mediated models of 
supersymmetry breaking; note that in this case the de- 
caying particle does not dominate the energy density of 
the Universe. Gravitinos with mass ^ 100 GeV — 3 
TeV have a lifetime r > 1 s, which can ruin the success of 
nucleosynthesis if gravitinos are produced in large abun- 
dances 1^. This, as well known, leads to constraints on 
the reheating temperature of the Universe Tr ||]. On 
the other hand, for TO3/2 > 10 TeV the gravitino decays 
before nucleosynthesis and has no effect on the light ele- 
ment abundances^. However, even in this case a signifi- 
cant upper limit on Tr can be derived from the following 
argument. If the LSP has a mass > 100 GeV, gravitinos 
(as heavy as 10^ GeV) decay much after the freeze-out 
of LSP annihilation. The overclosure bound then results 
in a constraint on the gravitino number to entropy ratio, 
< 4x 10""'^^, even when the gravitinos decay before 
the onset of nucleosynthesis. For thermal gravitino pro- 
duction, where rig/s/s ~ 10-"(Tr/10" GeV) §, this 
results in the limit Tr < 10^^ GeV; including possible 
non-thermal production of gravitinos |^] will presumably 
sharpen this limit. 

Finally, let us point out possible implications of deviat- 
ing from two major assumptions which we made through- 



gravitino mass which is much larger than superparticle 
masses can be naturally found, for example, in models of no- 
scale supergravity [pO{. 
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out this article: a significant branching ratio of primary 
infLaton decays into known SM particles (possibly includ- 
ing their superpartners), and allowed (pair-)production 
of X particles in the scattering of matter particles. First, 
consider the case when the inflaton exclusively decays 
to exotic light particles while x is produced through its 
coupling to matter particles. Assume that these exotic 
particles only couple with strength a' <C a to SM par- 
ticles. This results in a smaller Tmax! depending on the 
details of the model, one has to replace the factor 
in eq.(|^) either by aa'"^ or by a'^. We consider the lat- 
ter possibility and assume for simplicity that Tr remains 
unchanged. The thermal contribution will be signif- 
icantly reduced if the new value of Tmax is below by 
more than one order of magnitude or so. We therefore 



require that 



Let us find Vl'^''^^^ by assuming 



that Vl^-^ ~ 1. The main difference is that x production 
from inflaton decay now occurs through six body final 
state diagrams, with in eq.(^l]) replaced by a'^ . The 
requirement that Tmax > "^x results in a bound on a'^, 
see eq.(||). For rriy^ = 2 • IO^cx^'^Tyi, which saturates il^'^, 
and after using the bound on a'^ in (Eq), we find 



^dccay > 5 . -^q3 ^,32/21 



5/3 



m^/^MpiancklGeV 



(44) 



for X production in six body decay of the inflaton. Even 
for the most conservative choice, ~ 10"^^ GeV, eq.(p^) 
requires < 7 ■ 10^ GeV for ~ 0.01. Note that 
X production from hard-soft and hard-hard scatterings 
also only occurs through higher order processes, thus fl^ 
and are suppressed by a few orders of magnitude 
compared to expressions (p^), ( p^ and (p5|). Moreover, 
Tmax ^ TTi^ if a' is very small, implying that Q'^^ will be 
exponentially suppressed. In this case x production from 
inflaton decay will dominate, since it is only suppressed 
by a factor a'^. 

Second, consider the case where the inflaton decays to 
matter but x particles can only be produced from scatter- 
ings of some intermediate particle x' whose interactions 
with matter has a strength a' <C a. If a' is not very 
small x' will be in thermal equilibrium with matter and, 
provided Tmax > itt-x i ^^^^ ^^^Y be reduced by a statis- 
tics factor, since now only a small fraction of all soft-soft 
scatters can produce x particles. On the other hand, 
and will decrease much more. The hard-soft and 
hard-hard scatterings now can produce x's only in four- 
body final states, e.g. // x'x'xXi with cross section 
oc a^a'^. By using the expression in (^) a bound on a' 
is found in order for x' to be at thermal equilibrium for 
T > rriy-, which requires a'^m^ > H. If we now assume 
~ 1, taking g^, — 200, we have 



^docay > -^Q„^ 



m 



m 



X 



5„.12/7 ""X 

^ TO^Mpianck 1 GeV ' 



(45) 



from six body decay of the inflaton. Again taking 
= 0.01 and = 10^"^ GeV, this results in the bound 



TTT-x — 10^° GeV in order not to overdose the Universe. 
As a closing remark, we shall notice that such a rather 
contrived scenario might be realized for exotic x parti- 
cles (e.g. "wimpzillas"), but not when x is the LSP or 
a charged stable particle, which have electroweak gauge 
couplings to ordinary matter. 



IV. SUMMARY AND CONCLUSIONS 

In this article we studied the thermalization of pertur- 
bative inflaton decay products, with emphasis on appli- 
cations to the production of massive stable or long-lived 
particles x- We found that a thermal plasma should form 
well before inflaton decay is complete, if the theory con- 
tains light or massless gauge bosons; 2^3 reactions 
where a fairly energetic, nearly coUinear particle is emit- 
ted in the initial or final state play a crucial role here. 
The existence of light gauge bosons is required since only 
gauge boson exchange in the t— or m— channel leads to 
cross sections that significantly exceed a'^ /rn^, where a is 
the relevant (gauge) coupling strength and the infla- 
ton mass. This indicates that reheating might be delayed 
greatly if some scalar field(s) break all gauge symmetries 
during this epoch, which may naturally happen in the 
presence of supersymmetric flat directions. 

Even if massless gauge bosons exist, thermalization 
takes a finite amount of time. As a result, the maximal 
temperature of the thermal plasma will usually be well 
below (but can exceed the reheat temperature 
significantly); this limits the region of parameter space 
where thermal x production can play a role. On the 
other hand, it allows the very energetic primary infla- 
ton decay products to produce x particles either in colli- 
sions with the thermal plasma ("hard-soft" scattering), 
or with each other ("hard-hard" reactions). We esti- 
mated the rate for these reactions in the simple approx- 
imation where the primary inflaton decay products have 
energy for one thermalization time, and then have 

energy T. If is rather close to this will probably 
overestimate the true rate, since then the energy of the 
primary decay products will drop below the production 
threshold faster than in our approximation. On the other 
hand, if 3> m^, our approximation will likely be an 
underestimate, since in the process of thermalization a 
single particle with E ~ can produce several (most 

likely nearly coUinear) particles with > E > m^, 

all of which can contribute to nonthermal x production; 
note that this also allows x production from the scatter- 
ing of "hard" particles even if the primary inflaton decay 
products do not couple directly to x- We found that the 
hard-soft contribution will dominate over the hard-hard 
one if it is still kinematically allowed at T = Tr, but 
can otherwise be subdominant; either of these two new 
production mechanisms can exceed the rate from purely 
thermal x production. 

We also discussed x production in inflaton decay. We 
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pointed out that decays of this kind must be allowed, at 
least in four-body final states, if x particles can be pro- 
duced in collisions of primary inflaton decay products; in 
certain (somewhat contrived) scenarios one may have to 
consider six-body final states. In fact, in most cases this 
seems to be the most important nonthermal (but per- 
turbative) mechanism producing massive particles with 
m-^ < this contribution also often exceeds ther- 

mal X production by several orders of magnitude, even 
if is below the maximal temperature of the thermal 
plasma. Nonthermal x production therefore significantly 
sharpens limits on model parameters that follow from up- 
per bounds on the x relic density. For example, if is 
well above visible sector superparticle masses, each infla- 
ton decay will produce 0(1) lightest superparticle (LSP). 
If these LSPs were not in equilibrium at Tr , the bound 
^^LSP < 1 then implies that the reheat temperature must 
be at least ten orders of magnitude below the inflaton 
mass, i.e. the inflaton decay width must be at least 26 or- 
ders of magnitude smaller than m^. This does not seem 
to be very plausible. As well known, the requirement 
^^LSP < 1 imposes severe constraints on the parameter 
space of many supersymmetric models if the LSP was in 
equilibrium at Tr [ pl| . Our analysis indicates that it is 
quite difficult to evade these constraints by changing the 
cosmology. Similarly, we found that stable charged par- 
ticles can only be tolerated if they are too heavy to be 
produced in inflaton decays. 

Many of the results presented in this paper are only 
semi-quantitative. Unfortunately in most cases signif- 
icant improvements can only be made at great effort. 
For example, a more accurate treatment of thermaliza- 
tion would require a solution of the Boltzmann equations 
in the presence of a non-trivial, but non-thermal, back- 
ground of relativistic particles. Once a thermal plasma 
has been established, a proper treatment of the slow- 
down of primary inflaton decay products would require 
a careful treatment of the full momentum dependence of 
the particle distribution functions. On the other hand, 
our estimates of inflaton decay branching ratios should 
be quite reliable if rrix > Tr, (which is required for x 
not to have been in thermal equilibrium at Tr); even for 
many-body decays, details of the matrix elements should 
change our estimates only be 0(1) factors. Fortunately 
we found that this is often also the most important of 
the new, non-thermal mechanisms for the production of 
massive particles at the end of inflation. We therefore 
conclude that cosmological constraints on models with 
stable or long-lived massive particles are (much) more 
severe than had previously been thought. 
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